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Abstract 

It is shown how the quantumness of atom-atom correlations in a trapped bosonic gas can be made ob- 
servable. Application of continuous feedback control of the center of mass of the atomic cloud is shown 
to generate oscillations of the spatial extension of the cloud, whose amplitude can be directly used as a 
characterization of atom-atom correlations. Feedback parameters can be chosen such that the violation of a 
Schwarz inequality for atom-atom correlations can be tested at noise levels much higher than the standard 
quantum limit. 
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In optics correlations between photons in a light field have been demonstrated using a Hanbury 
Brown Twiss setup 111. There the incoming light beam is split by a semitransparent mirror and 
the two outputs are measured with photodetectors. By displacing detectors normal to the detected 
beams and by introducing a variable time delay, from the coincidences of detections the intensity- 
intensity correlation of the quantized light field, 

C(1,2) = (S/(1)/(2)S). (1) 

is obtained, where 1 and 2 represent the space-time coordinates (xi, ti) and (x2, ^2), respectively. 
Here / = E^~'^E^^^ with E^'^'^ being positive- and negative-frequency parts of the (here scalar) 
electric field of the light beam and o o denotes time and normal-ordering. 

Interpreting the expectation value in Eq. dU) as being based on a proper probability density, 
application of the Cauchy-Schwarz inequality leads to 

C(l,2)< v'C(l,l)C(2,2). (2) 

This inequality is based on the assumption of a classical random process with correlations de- 
scribing either complete randomness or bunching of photons fl. However, quantum states of 
light exist that cannot be described in this way and violate the inequality @. This has been used 
as a criterion for true quantumness — or non-classicality — of the incoming light beam, and has 
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been experimentally observed as antibunched light |yl|4llSla]- 

These concepts can also be applied to ultracold bosonic gases, where the electric field is re- 
placed by the matter- field (f){x,t). However, different from optical correlations measured by 
absorbing photodetectors, atom correlations will not occur in normal operator ordering. Never- 
theless, an analogy between photonic and atomic correlations can be drawn. In fact, atom-atom 
correlations have been recently measured in a way analogous to the Hanbury Brown-Twiss optical 
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setup ILZHSIISD- Furthermore, there have been also approaches where the measurement of losses 
due to three-body recombinations has been used to infer on those correlations 

In these techniques either the gas has to be released from the trap or atomic losses play a major 
role. One might, however, manipulate or drive the system in such a way as to map the internal 
atom-atom correlations into easily accessible mesoscopic observables, that may be observed in 
situ without loosing atoms from the trap. This is the approach taken in this Letter: By applying a 
continuous feedback control of the atomic cloud's center of mass and by observing the cloud size, 
it is shown that atom-atom correlations can be detected and used to test for the quantumness of 
atom-atom correlations, in close analogy to antibunching in the case of photons. 



Consider an ideal bosonic gas with atoms of mass m, that is kept in a harmonic trap of 
frequency u. A feedback loop is continuously applied to compensate for the motion of the center 
of mass of the gas X, i.e. to damp the collective motion of atoms. It consists of the continuously 
repeated application of a measurement of X with measurement outcome X„^ and resolution ctq, 
and a corresponding shift by — ^o-^m- Both processes shall be much faster than the free oscillation 
of the system with trap frequency u, so that an instantaneous action can be assumed. 

Experimental implementations of such a feedback loop are based on the collective interaction 
of the atoms with far off-resonant optical probe fields. Two-photon transitions, whose strength 
depend on the positions of atoms in the probe field, then lead to a redistribution of intensities in 
optical-field modes. The latter may be detected and used as an input signal for a subsequent control 
action using optical phase shifts of the same probe field. Such a scheme has been realized with 
atoms in optical lattices flil3i1 or may be realized as an extension of the single-atom experiment 
by Fischer et al. Possibly it may also be integrated into the magneto-optical or optical trap 
configuration, by detecting and modulating the trapping laser fields. 

Given the continuous application of feedback at rate 7, two parameters determine the feedback 
dynamics: the rms time-integrated measurement resolution a = ao/ ^ and the feedback shift rate 
C = (^o7- The time evolution of the A^-atom density o perator of the system, ^at, is then described 
by the master equation of quantum Brownian motion 
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dt QN = -^[H, Qn] + {X, qn}] - ^[X, [X, ^;v]] - ^[A [P, Qn]]- (3) 



Here H is the Hamiltonian of the non-interacting atomic gas in the trap potential and the center- 
of-mass (cm) operator and its canonically conjugate total momentum read 

X = ^ J dx (f)\x) X (j){x) , P = —ih J dx 4>\x) dx (i){x) , (4) 

with the atomic field (j){x) obeying the bosonic commutator relation [(j){x) , (j)^ {x')] =5{x—x'). 

The stationary behavior of the cm, obtained from Eq. is an exponential damping at rate C/2 
of the coherent oscillation: \imt-,oo{X{t)) = 0. Moreover, its rms spread converges exponentially 
at the same rate to the non-zero stationary value 
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lim V([AX(t)]2) = AX„ (5) 



where AX = X—{X). It represents the noise left in the cm after a time of the order of C ^ needed 
for damping the coherent cm oscillation. This noise is determined solely by the parameters of the 



feedback and trap: 

AXs = 6X0 v/(r/+r/-i)/2, (6) 
with the rms spread of the cm in the ground state of the trapping potential being 



6X0 = ^h/{2Nmuj), (7) 

which serves as the standard quantum limit (SQL) for the cm coordinate. The parameter 

rj = 6X0^ {(a') (8) 

specifies the ratio of spatial localization due to the potential over that due to the feedback. Note, 
that Eq. Q attains the SQL as a minimum value for 77 = 1 but is otherwise much larger. 

Our goal is to describe atomic correlation effects in the dynamics of the atomic density. For 
that purpose we need the single- atom density matrix 

p{x,x',t) = {^^{x')^{x))t, (9) 

with {. . .)t = Tr[. . . g^it)]- To obtain the dynamical evolution of this density matrix from Eq. ©of 
course is prevented by the correlations in the many-atom systems. That correlations play a role in 
Eq. dSj) can be seen from the occurrence of products of operators X and P, that contain products 
of four field operators [cf. Eq. ©J, similar to atom-atom interactions. However, recently it has 



been shown 111 811 . that despite of these problems, the single-atom density matrix can be obtained 
via a procedure that we may briefly outline here: 

Instead of the single-atom density matrix, the joint Wigner function of single atom (variables 
X, p) and cm of the other A^ — 1 atoms (variables X, P) is considered: 

W{x,p-X,P,t) = {27rhy^ J dx'e-''''^/^ J dX' J dP' 
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From Eq. © a closed Fokker-Planck equation follows for this distribution, which is of linear type 
with positive semi-definite diffusion matrix, leading thus to a bound analytic Green function of 
Gaussian type lioi . In consequence, given the initial conditions, analytic solutions for this Wigner 
function can be obtained, from which the solution for the single-atom density matrix are derived 
by integration over the auxiliary phase- space variables: 

p{x+x',x-x',t) = JdxJdPjdpW{x,p;X,P,t)e'^'P'''. (11) 



Thus in principle the complete atomic density profile could be obtained. Here we focus on the 
rms spread of the corresponding atomic density, 



^x{t) = I J^x"^ p{x,x,t) - J^xp{x,x,t) 
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(12) 



giving us information on the quantum-statistically averaged temporal evolution of the extension of 
the atomic cloud. From a complete solution of the Fokker-Planck equation for dTOb this variance 
can be shown to exponentially converge at rate C/2 to the asymptotic behavior limt^oo Ax(t) ~ 
Axa(t), defined by 



Axa(t) = VAXs^ + a,2(t). (13) 

In this equation the first term in the square root is given by the constant stationary rms cm 
spread [cf. Eq. ®], whereas the second term is explicitely time dependent and reads 

^.'W = W^)]mx))^- J^{4>\x)q{x,t)^{x)^\x')q{x',t)^{x'))^, 

(14) 

with the expectation value being taken with respect to the initial A^-atom density operator ^Ar(O). 
The explicit time dependence of Eq. (fT4h is given by the single-atom quadrature, defined as 

ifiO 

q(x,t) = xcos(ujt) -sm(ujt). (15) 

mu 

Equation (fT4b represents the central result of our Letter. In its second part it contains an atom- 
atom correlation function with four matter- field operators, showing that due to the feedback these 
correlations show up in the observable mesoscopic size of the cloud, cf. Eq. dTSt . Thus the atom- 
atom correlations will become visible in a purely single-atom property. 

Due to its explicit time dependence, in general there will be no stationary size of the atomic 
cloud, but instead the cloud will periodically breath. It should be emphasized, that this breathing 
has nothing in common with the well-known collective oscillations of a Bose gas or condensate. 
The latter rely on the presence of atomic collisions, whereas the breathing discussed here is an 
effect solely produced by the feedback. Feedback of course effectively mediates interactions be- 
tween atoms, so that a certain analogy to true atom-atom interactions can be drawn. However, 
whereas collisions of ultracold atoms lead to Hamiltonian terms, the feedback in addition provides 
non-unitary parts of the time evolution [cf. Eq. ©]. 

Thus, after a transient behavior during a time of the order of the cloud size will oscillate at 



twice the trap frequency, which resembles single-atom quadrature squeezing ll20ll . The amplitude 



of this oscillation will of course depend on the initial quantum state and its correlations at time 
t = 0, before the feedback has been turned on. Since the SQL for a single atom in the trap reads 
6x0 = 6Xq^/N, it follows that the atomic cloud size (fT3l) can become smaller than that value. The 
condition for quadrature squeezing (QS) on the single-atom level would then be that at some time 
during the half period n/u the following inequality is fulfilled: 

Axs,{t) < 6x0 (single-atom QS). (16) 

However, there is more to Eq. (O than single-atom QS. It can be revealed by taking a closer 
look at the structure of Eq. (fT4l) and applying the ideas developed in the context of photon anti- 
bunching. The probability density for finding an atom at position x is undoubtedly defined as 

P{x) = {4>\x)4){x))/N. (17) 

A quasi joint probability density for two atoms being at positions x and x', that is consistent with 
the definition dlTt . reads 

P{x,x') = {4)\x)4){x)(j)\x')(j){x'))/N\ (18) 

Note that different from the optical case this correlation is not normally ordered. Consistency 
means here that for N atoms the marginals of P{x, x') reproduce the correct probability density: 

Jdx' P{x,x') = P{x), JdxP{x,x') = P{x'). (19) 

Clearly P(x,x') is not a proper probability density in general. However, when interpreting the 
atomic fields as classical ones: — > it becomes a proper classical joint probability density. 
In this classical interpretation, we may now apply the Schwarz inequality that states 

Jdxjdx' q{x)q{x')P{x,x') < jdxq^{x)P{x) (classically). (20) 

Applying this result to Eq. (fT?l) one obtains a classical inequality for the contribution to the size 
of the atomic cloud: o'^it) > 0. In consequence, via Eq. dlSt one arrives at the classical inequality 
for the size of the atomic cloud: Axa(t) > AXg. 

A violation of the Schwarz inequality would indicate true quantum correlations between atoms 
as opposed to classical ones, since then the second expectation value in Eq. (fT4ll cannot be de- 
scribed by integration over a proper probability density. The condition for this case is then 

Ax^{t) < AXs (Schwarz violation). (21) 
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Indeed, since Eq. (fT4l) contains atom-atom correlations in the second term on the rhs, one may 
interpret a violation of the Cauchy-Schwarz inequality, as formulated in Eq. d^TT) . as a test for true 
quantum correlations between atomic pairs. Note, that the concept of defining quantumness here is 
the same as in the case of photon antibunching in the context of the optical Hanbury Brown-Twiss 
experiment. 

A Schwarz violation at time t can be easily shown to be equivalent to a violation of the (classi- 
cal) inequality 



where Aq and AQ are the rms spread of the "single-atom" quadrature, corresponding to the atomic 
density, and the rms spread of the cm quadrature, respectively. The nature of such a violation can 
be understood by considering the special case where the Schwarz violation occurs at a time where 
the quadrature ([T51 reduces to the atomic position, i.e. q{t) ^x. The relation (E^ then states, that 
classically the size of the atomic cloud Ax is always equal or larger than the rms spread of the cm 
of the cloud. In other words, the cm of the object is well localized within the spatial extension of 
the cloud. A violation of this classical inequality would then correspond to cases where the cm 
coordinate reveals a rms spread larger than the cloud's size. In the extreme case, the atomic cloud 
can then be seen as an almost pointlike object, its internal distribution not being resolved, whose 
(cm) coordinate fluctuates. 

It is thus the transition from an atomic cloud with well localized cm to a quasi pointlike object 
with large fluctuation of its coordinate, that corresponds to a transition of classical to quantum 
atom-atom correlations in the cloud. Clearly for the general case, i.e. a Schwarz violation at an 
arbitary time withing the half period n/uj, the above interpretation correspondingly holds for a 
specific quadrature q{t) instead of position. 

In overall thus two levels of quantum or - if one wishes to use this term - non-classical behavior 
can be distinguished by the amplitude of the oscillation of the atomic-cloud size. These two 
however, single-atom QS and Schwarz violation, do not form a unique hierarchy: For different 
parameters i] of the feedback mechanism the two boundaries [cf. Eqs (fT6t and (12111 1. appear in 
different orders. The order depends on whether AX^ < 5xo or not. The parameter range for this 
condition is obtained as 



In the case ^ oo this range includes all possible values of r] and thus for a truly macroscopic 




(22) 



AX, < 6x0 for N-VN^-1 <V< N+VN^-l. 



(23) 
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FIG. 1: Observable size of the cloud versus time: For the parameter range given by Eq. ( 123 1 (left part) the 
Schwarz boundary is below the one for QS whereas for values of t] outside this range (right part) the reverse 
is true. 

system, before a Schwarz violation can be observed always first single-atom QS appears. This 
case is depicted in the left part of Fig.[T] 

However, for a finite system values for 77 can be found outside the range (l23t which then re- 
veals an exchange of the order of single-atom QS and Schwarz violation. In this case quantum 
correlations between atoms can be observed without the presence of single-atom QS, cf. right part 
of Fig.lH The corresponding parameter range is given by 



In this range of 77 apparently the quantumness of correlations can be detected at much higher noise 
levels than the single-atom SQL Sxq, which may render its experimental observation substantially 
more feasible. Moreover, the two ranges in Eq. (l24l) correspond to weak and strong feedback 
localization, respectively, and thus may allow for a suitable combination of values for the shift rate 
C and the rms time-integrated measurement resolution a in an experiment. Last but not least this 
scheme may be applied even for the extreme case of only two indistinguished bosonic atoms, for 
which the range of possible values of 77, according to the case (l24ll . becomes even more broad. 

The criterium for quantumness of correlations can be tested directly from measurements of 
the size of the atomic cloud Ax{t) over half a period of the trap oscillation, after the system has 
reached its asymptotic behavior within a delay time of the order of l/C- Experimentally a suffi- 



r] < N-VN^-1 or rj > N+VN^-1. 



(24) 
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ciently large number of sequences of feedback evolutions of varying time duration and final cloud- 
size measurements have to be performed. Each sequence starts with the identically prepared ini- 
tial quantum state, whose correlations are to be detected. Thus the final cloud-size measurements 
can be arbitrarily destructive and can be performed for example by density-profile measurements. 
These may be implemented by absorption imaging lllj . dispersive light scattering |22], or pos- 
sibly phase-modulation spectroscopy 11231] . Thus together with possible experimental techniques 
to generate the required feedback loop lisifl^ . an experimental implementation of the presented 
scheme and a test for quantum atom-atom correlations in bosonic gases seem to be feasible. 

In summary we have shown that the quantumness of atom-atom correlations in a trapped 
bosonic gas can be made observable as size oscillations of the atomic cloud via feedback. For 
weak and strong feedback localization a Schwarz violation for atom-atom correlations of a gas 
with finite atom number can be observed in the absence of single-atom QS at correspondingly 
higher noise levels than the SQL. Together with the feasibility of implementing this scheme with 
present experimental techniques, this may allow for detecting the quantumness of atom-atom cor- 
relations in a bosonic gas. 
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